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Abstract 

Polyelectrolyte gels are a very attractive class of actuation materials with remarkable elec¬ 
tronic and mechanical properties with a great similarity to biological contractile tissues. 
They consist of a polymer network with ionizable groups and a liquid phase with mobile 
ions. Absorption and delivery of solvent lead to a large change of volume. This mechanism 
can be triggered by chemical (change of salt concentration or pH of solution surrounding 
the gel), electrical, thermal or optical stimuli. Due to this capability, these gels can be used 
as actuators for technical applications, where large swelling and shrinkage is desired. In the 
present work chemically stimulated polymer gels in a solution bath are investigated. To ad¬ 
equately describe the different complicated phenomena occurring in these gels, they can be 
modeled on different scales. Therefore, models based on the statistical theory and porous 
media theory, as well as a coupled multi-field model and a discrete element formulation are 
derived and employed. 

A refinement of the different theories from global macroscopic to microscopic are pre¬ 
sented in this paper: The statistical theory is a macroscopic theory capable of describing 
the global swelling or bending, e.g., of a gel film, while fhe general fheory of porous me¬ 
dia (TPM) is a macroscopic confinuum fheory which is based on fhe fheory of mixtures 
extended by the concept of volume fractions. The TPM is a homogenized model, i.e., all 
geometrical and physical quantities can be seen as statistical averages of the real quantities. 
The presented chemo-electro-mechanical multi-field formulalion is a mesoscopic fheory. 
If is capable of giving fhe concenfrafions and fhe elecfric pofenfial in fhe whole domain. 
Finally fhe (micromechanical) discrefe elemenf (DE) fheory is employed. In Ibis case, fhe 
confinuum is represenfed by disfribufed parficles wifh local inferacfion relafions combined 
wifh balance equafions for fhe chemical field. This mefhod is predesfined for problems 
involving large displacemenfs and sfrains as well as disconfinuifies. 

In Ibis paper, fhe coupled mulli-field model and fhe discrefe elemenf model for chemical 
sfimulafion of a polymer gel film wifh and wifhouf domain deformation are employed. 
Based on fhese resulfs, fhe presented formulalions are compared and conclusions on fheir 
applicabilify in engineering pracfice are finally drawn. 
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1 Introduction 

Polymer gels are ionie eleetroaetive polymers. They eonsist of a polymer network 
with eharged groups and a liquid phase with mobile ions, see Figure 1. By a ehange 
of the ehemieal milieu in the gel surrounding solution or by an applied eleetrie 
field, the swelling of these gels ean be triggered. Due to this eapability, eleetrolyte 
polymer gels are designated as actuators for technical applications where large 
swelling and shrinkage is desired, such as for artificial muscles or other chemo- 
electro-mechanical actuators. The swelling behavior of these gels results from the 
equilibrium of different forces such as osmotic pressure forces, electrostatic forces 
and viscoelastic restoring forces. 
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electrolyte polymer gel solution 

Figure 1. Polymer gel network with bound anionic groups and mobile ions in solution 


2 Modeling of Polymer Gels 

To describe the different phenomena adequately, the modeling is conducted on dif¬ 
ferent scales. Starting on an atomistic level, the crystal structure of the material can 
be resolved; by increasing the length scale, the structure can be modeled by molec¬ 
ular dynamics or particle models which take into account the effects occurring on 
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the microscale. By further increasing the scale, the whole structure can be investi¬ 
gated on the meso- or macroscale, e.g., when the prediction of the global failure of 
the structure is desired. 

In the following, some key points on the statistical theory, the porous media theory, 
the coupled multi-field model and the discrete element theory are given. 


2 .1 Statistical Theory 


The statistical theory is a method to describe the global swelling ratio of polyelec¬ 
trolyte gels. This macroscopic theory is based on the swelling of network structures 
by Flory [1]. The conditions in the gel and of the gel surrounding solution is de¬ 
scribed by a change of the Gibbs free energy AF. The total free energy is the sum 
of the free energies of mixing AFm , of elastic deformation AFei and of the mo¬ 
bile ion phases AFion ion inside and outside the gel, see [2, 3, 4]. The equilibrium 
state is characterized by a minimum of the total free energy difference AF. In the 
standard form, this theory is valid for chemical and also thermal stimulation. More 
details on this theory together with numerical simulations can be found in the ref¬ 
erences given above or, e.g., in Wallmersperger et al. [5, 6] or Maurer et al. [7, 
8 ]. 


2.2 Porous Media Theory 


The general theory of porous media (TPM) is a macroscopic/mesoscopic contin¬ 
uum theory which is based on the theory of mixtures extended by the concept of 
volume fractions. In this theory neither the local porous micro structure nor the 
actual geometrical distribution of all the constituents have to be known [9, 10]. 

The TPM is formulated by conservation equations for the different phases or mix¬ 
tures. It comprises local mechanical and chemical unknowns. The TPM is applica¬ 
ble first of all for the chemical stimulation but can be extended easily to other kinds 
of stimulation by adding new equations and coupling terms. More details can be 
found, e.g., in Ehlers et al. [11] and Kunz [12]. 


2.3 Chemo-Electro-Mechanical Multi-Field Formulation 


The coupled multi-field formulation [13, 14] is formulated by different balance 
equations. It is capable of giving the local chemical, electrical and mechanical un¬ 
knowns. 
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The chemical field is described by a convection(-migration)-diffusion equation for 
the different mobile species. The electric field is formulated by the Poisson equa¬ 
tion. In contrast to the TPM, the electric field is directly obtained by solving the 
Poisson equation, i.e., no additional conditions prescribing the jump between gel 
and solution have to be used. 


The mechanical field is formulated by the momentum equation. Normally, the de¬ 
formation process is very slow, i.e., the influence of the inertia term is very small 
and the second order in time contributive term can thus be neglected. 


The chemical and electrical field are solved simultaneously. The different concen¬ 
trations in gel and solution give rise to an osmotic pressure term in the momentum 
equation leading to a change of the mechanical displacement. 
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Figure 2. Modeling strategies for electrolyte polymer gels 


2.4 Discrete Element Theory 


The micromechanical behavior is described by the discrete element (DE) theory. 
In this model, the continuum is represented by distributed particles comprising a 
certain amount of mass. The particles interact with each other mechanically [15] 
by a superimposed truss or beam network of massless elements [16]. 

The mechanical behavior, i.e., the dynamics of the system, is obtained by solving 
the Newton’s equations of motion 


mxi = ^ Fij 
j 

Qui = Mij 
j 


( 1 ) 

( 2 ) 
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where x is the aeeeleration of the partieles (in x-direetion), Cj the angular aeeeler- 
ation, m the mass of the partiele, and 0 the rotational moment of inertia. Fij and 
Mij are the generalized forees and moments, respeetively. 

Additionally to the meehanieal field, the ehemieal field, i.e., the ion movement in¬ 
side the gel and from the gel to the solution, is deseribed by diffusion equations for 
the different mobile partieles. The eoneentration differenees between the different 
regions of the gel and between gel and solution form the osmotie pressure differ- 
enee, whieh is a main eause for the meehanieal deformation of the polyeleetrolyte 
gel film. 


3 Numerical Simulation 


In this seetion, the ehemieal stimulation, i.e., the ehange of salt eoneentration, of 
polymer gels in a solution bath is investigated. 

In the first test ease, the ehemo-eleetrieal behavior of a polyeleetrolyte gel film 
(5x5 em^), fixed in a solution bath is eonsidered (see Figure 3), while in the seeond 
test ease the domain deformation is taken into aeeount. In both 2D simulations, 
the eoneentration of the bound anionie groups in the gel is fixed. At f = 0 the 
eoneentration in the solution has been suddenly deereased from Cs,o = 2mM to 
Cs = ImM. The parameters used in the numerical simulations are given in Table 1. 



3.1 Chemo-Electrical Simulation on a Fixed Domain 


In this test case, only the chemo-electrical field will be investigated without con¬ 
sidering any domain deformation. 

The ion fluxes in the whole fixed gel-solution domain and the resulting time-history 
of the mobile concentrations and the potential difference between gel and solution 
are computed by using the coupled chemo-electric multi-field formulation. 

With progressing time, the concentrations of Na^ and Cl~ in the gel are decreasing 
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Domain length: 

lx = 15 mm 

Init. coneentration of 

anionie 

c^- = 2mM 



rest groups: 


,0 

Domain width : 

ly = 15 mm 

Coneentrations in the solution: 

= Cs = ImM 

Gel length: 

Ixg = 5 mm 



= Cs = ImM 

Gel width: 

lyg = 5 mm 

Initial eoneentrations in 

lution: 

the so- 

Cs,o = 2mM 

Temperature: 

T = 293 K 

Diffusion eoefficients: 


D = 10-^m2/s 


Table 1 

Parameters for the numerieal simulation 

from the initial values n = 3.2361mM, = 1.2361mM) until the steady 

state (c^a+ ~ 2.4142mM, = 1.4142mM) is reaehed. In Figure 4 the station¬ 

ary solution of the ion eoneentrations and the eleetrie potential in the gel and the 
solution is depieted. The eleetroneutrality eondition is fulfilled in both areas (gel 
and solution) outside the small boundary layer region. 


3.2 Discrete Element Simulation 


In the seeond test ease, the swelling behavior of a polymer gel fiber in a solution 
bath is investigated with the diserete element (DE) method; i.e., the eoupled ehemo- 
eleetro-meehanieal field is solved under eonsideration of the domain deformation. 

fs) 

For t < 0, the anionic gel (c^- = 2mM) is located in a solution bath (c^ = = 

= 2mM) in which the stationary equilibrium is fulfilled. At t = 0, the concen¬ 
tration in the solution bath is reduced to Cs = ImM. 

In the DE simulation, only the gel - i.e., without surrounding solution - is investi¬ 
gated. The boundary conditions at the gel/solution interface are computed at each 
time step by solving the neutrality condition in gel and solution 

E(C«C'a) = 0 (3) 

Q. 


where Ca is the charge of the ion species a - together with the Donnan-equilibrium 


^Na+ 

^NA+ 


(s) 

^CL- 

M 

'^ci- 




39) 


Sq) 

^ci- 


(4) 


resulting from identical electro-chemical potential in both regions. The resulting 
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Figure 4. 2D-steady state results of the ion concentrations {Cl~ (bottom left), Na^ (bottom 
right)) and the electric potential if; in the gel and the solution without applied external 
electric field. 

osmotic pressure Att 

Nf 

(5) 

CK=1 


due to the eoneentration differenees leads to a ehange of the swelling ratio and 
therefore of the gel geometry. R is the universal gas eonstant, T the temperature 
and and are the eoneentrations of the speeies a in the solution and the gel, 
respeetively. 


Due to the fixed number of anionie eharges, the number of moles in the gel remains 
identieal, i.e., the eoneentration of bound groups varies by 



( 6 ) 


where V is the aetual volume of the gel and Vq the volume in referenee state. 

Sinee the boundary eonditions are dependent on the salt eoneentration in the solu¬ 
tion Cg and the aetual bound eoneentrations in the gel c^?, different values for the 
eoneentrations at the boundaries have to be preseribed for eaeh time step. 

In Figure 5, the initial eonditions of the bound eharges (top left) and of the mo- 
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bile anions (top right) are depieted. Due to the reduetion of the mobile ions in the 
solution, a new boundary condition at the gel/solution interface is prescribed and 
thus, the concentration of mobile anions in the gel is also decreasing. A typical plot 
of the mobile anions in the gel is given at a specific time see Figure 5 (middle 
right). Due to the resulting swelling of the gel film, the bound anionic concentra¬ 
tion is reduced, see Figure 5 (middle left). The stationary solution of the bound 
anionic charges and mobile ions versus the new gel geometry is depicted in Figure 
5 (bottom). Note, the cationic concentrations in the gel may be calculated by 


c 


(ff) 

Na+ 



+ Ca- 


(7) 
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Figure 5. Fixed negative charges (left) and mobile negative ions (right) versus the gel ge¬ 
ometry at i < 0 (top), t = ti (middle) and t —)• oo (bottom). 


The swelling of the gel is depicted in Figure 6: the initial geometry is represented 
by the black line while the swollen gel in steady-state is given by the grey one. 
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Figure 6. Swelling of the gel fiber: gel geometry at f < leqO (black) and steady-state 
solution (grey). 

3.3 Comparison 


In order to compare both test cases, in Figure 7, the concentrations of the fixed 
charges and the mobile ions in steady state without domain deformation (bottom 
left) and with domain deformation (bottom right) are depicted. The initial concen¬ 
trations at f < 0 versus x are given in Figure 7 (top). 

It can be seen that due to the resulting enlargement of the domain (Figure 7 (right)), 
the bound concentration ca- decreases according to Eq. (6), while the total num¬ 
ber of moles of bound charges remains constant. Due to the decrease of the bound 
anionic concentration, new concentration boundary conditions are prescribed - see 
Eqs. (4) and (7) - and therefore, the concentration of mobile cations in the gel also 
decreases. Summarizing, the positive domain deformation decreases the concentra¬ 
tion differences; therefore the differential osmotic pressure decreases. This means, 
the enlargement of the gel domain is reduced by a kind of restoring force. 


4 Conclusion and Outlook 


Substantially different modeling strategies for the swelling behavior of polyelec¬ 
trolyte gels under electro-chemical stimulation were presented in this paper. 

If the global macroscopic, mechanic response is the only property of interest, the 
use of a simple statistical theory leads to sufficient results. Both the porous media 
theory and the multi-field formulation work within the assumptions and limitations 
of standard discretization schemes such as finite element or finite difference meth¬ 
ods. Due to the field formulations, continuous fields of mechanical unknowns like 
displacements, strains and stresses, and of chemical and thermal unknowns are ob- 
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Figure 7. Concentrations of the bound anions Cfix = c^- and the mobile cations 
c'^ = cjva+ anions c~ = Cci- versus the Xg -coordinate: initial state (top), steady 
state without swelling (bottom left) and steady state with swelling (bottom right). 


tained in a straightforward way. These formulations are valid for quasi statie and 
transient problems. 


In order to model the chemical stimulation, for the statistical, the porous media 
as well as the discrete element theory, the gel is investigated only; therefore the 
Donnan equilibrium at the boundary has to be prescribed and the electro neutrality 
condition in the whole domain has to be fulfilled. The coupled multi-field formula- 
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tion is applied for the whole gel-solution domain and the electroneutrality condition 
is not (directly) prescribed. 


In this paper the coupled multi-field model and the discrete element model for 
chemical stimulation of a polymer gel film with and without domain deformation 
have been employed and compared. It can be seen that the chemo-electrical multi- 
field formulation is an optimal method to give the chemial and electrical unknowns 
in the whole domain. If the chemo-mechanical excitation should be considered, 
the discrete element approach in its present form is also predestined. The advan¬ 
tages of discrete element representations are, e.g., the direct physical access to the 
system during all the time steps of the simulation, and the representation of large 
deformations and strains. 
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